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Abstract 

A first principle derivation is given of the neutrino damping rate in real- 
time thermal field theory. Starting from the discontinuity of the neutrino self 
energy at the two loop level, the damping rate can be expressed as integrals 
over space phase of amplitudes squared, weighted with statistical factors that 
account for the possibility of particle absorption or emission from the medium. 
Specific results for a background composed of neutrinos, leptons, protons and 
neutrons are given. Additionally, for the real part of the dispersion relation 
we discuss the relation between the results obtained from the thermal field 
theory, and those obtained by the thermal average of the forward scattering 
amplitude. 
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I. INTRODUCTION 



A neutrino propagating in matter do not longer respect the vacuum energy momentum 
relation. The modification of the neutrino dispersion relation is caused by their coherent 
interaction with the particles in the background and can be accounted in terms of an index 
of refraction |T| or an effective potential 0. The topic of neutrino propagation in matter 
became of prime relevance after Wolfenstein || study of the neutrino refractive index in 
matter, and later when Mikheyev-Smirnov j|] recognized the resonant neutrino flavor os- 
cillations triggered by matter effects. The Mikheyev-Smirnov effect has become the most 
popular explanation of the solar neutrino deficit ||. 

In general, the neutrino dispersion relation is a complex function uj k = oj(k). According 
to the Thermal Field Theory (TFT) matter contributions to the real and imaginary parts 
of u(k) arise from the temperature and density-dependent part of the neutrino self-energy. 
To leading order in (g 2 /M^) the real part of the dispersion relation is proportional to the 
particle-antiparticle asymmetry in the background. If the asymmetry is small or zero the 
imaginary part of uj(k) and corrections of order g 2 /M^ v to the real part may be important 
because they do not depend on the differences between the number of particles and antipar- 
ticles. This may be the case in the early Universe, when the medium was probably nearly 
CP- symmetric. 

Special attention has been given to the calculations of the 0(g 2 /M^r) corrections to 
the real part of the neutrino dispersion relations . Within the framework of the TFT 
these corrections arise from the momentum-dependent terms of the boson propagator in the 
self-energy diagrams. 

The imaginary part of the index of refraction for neutrinos propagating in a CP- 
symmetric plasma composed of electrons, neutrinos and their anti-particles has been consid- 
ered in references p|,P,[TDfl . These calculations have been based upon the computation of the 
neutrino reaction rates, assuming massless background fermions. In our opinion the relation 
of the neutrino damping rate to the self-energy discontinuties deserves further consideration. 
Additionally this work addresses the effects of the nucleons background contributions and 
the fermion mass correction. 

A systematic method to compute the damping rate from the imaginary part of the self 



energy can be formulated in terms of the Cuttosky thermal rules. Weldon [11] and Kobes 12 



proved that the imaginary part of the self energy can be organized in a form that includes 
square of amplitudes of the various processes obtained from the cuts of the self energy and 
weighted with the appropriated statistical factors. The examples discussed in those papers 
are always at the one loop level. As it shall be discussed, the calculation of the neutrino 
damping rate require to consider the self energy at the two loop level, the interpretation in 
terms of square of amplitudes of the allowed processes will be proved to remain valid. The 
approximations required to recover the results obtained from the optical theorem will be 
clearly stated. 

In this work we use the method of real-time thermal field theory to carry out a carefully 
calculation of the imaginary part of the neutrino dispersion relation in a medium composed 
of electrons, protons, neutrons, neutrinos and their anti-particles. As already mentioned 
the contributions to the imaginary part of the neutrino self-energy vanishes at the one loop 
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level, so we have to consider the contributions at the two loops level. 

The paper is organized as follows. In the next section we briefly review those ingredients 
of the FTF formalism that are required to accomplish our calculations. In Section |T| we 
re-derive the real part of the dispersion relation of a massless neutrino that propagate in 
a thermal background. We prove that utilizing the methods of finite -temperature field 
theory (at the g 2 order), the neutrino effective potential reduces to the thermal average 
of neutrino forward scattering amplitude. The calculation of the imaginary part of the 
dispersion relation is presented in Section |V[ The neutrino damping rate is extracted from 
the discontinuity of self-energy at the two loop level, it is expressed in terms of integrals 
over space phase of amplitudes squared, weighted with statistical factors that account for 
the possibility of particle absorption or emission from the medium. Specific results for 
a background composed of neutrinos, leptons, protons and neutrons are given. The last 
section contains a summary of our main results. 



II. BASIC FORMALISM 

The relevant quantity is the self-energy S, which embodies the effects of the background 
on a neutrino that propagates through it. According to the real-time formalism of the TFT 
p~3j^T5| , the real and imaginary part of S are given by the formulas 



Re £ = Re En, (2.1) 



2i( m (k-u) 1 -0(-k-u))' (2 - 2) 

where £ a b (a,b = 1,2) are the complex elements of the 2x2 self-energy matrix to be computed 
utilizing the Feynman rules of the theory. 9 is the step function and 

r]f(k ■ u) = [0(k ■ u)n f (x k ) + 9(-k ■ u)nf(-Xk)] , (2.3) 
where the thermal distribution is given by 

with Xk = /3(k ■ u — fif). Here, j3 is the inverse of the temperature and ///is the chemical 
potential associated with the fermion /. We have introduced the velocity four-vector of 
the background u^. In its own rest frame w M = (1, 0) and the components of the neutrino 
momentum k^ 1 are k^ 1 = (u, k). 

In the presence of a medium the chiral nature of the neutrino interactions implies that 



the self-energy of a massless (left-handed) neutrino is of the form [16| 

£(fc) = (a^ + bi)L, (2.5) 
where L — (1 — 7s)/2 and a, b are complex functions of the scalars 

lu — k ■ u, k = \J oj 1 — k 2 . (2-6) 
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In this case, the Dirac equation for the spinor U of the neutrino mode in the medium is 



E)U = 



(2.7) 



which has non-trivial solutions only for those values of u and k such that V 2 = 0, with 
= (1 — ajkp — bu^. Then, the dispersion relations uj k of the neutrino and antineutrino 



modes are obtained as the solutions of 



f{uJ K ,K) = 0, 



(2.8) 



and 



f(uj K ,n) = 0, 



(2.9) 



where 



f(uj,K) = (1 - a)(w -«)-&, 
f(uj, k) — (1 — a) (u + k) — b . 

In general, the solutions u K are complex, as usual we write 

uj k = u r - ij/2, 



(2.10) 



(2.11) 



where both u> r = Rea; K and 7 = —2 Imw K are real functions of n. A consistent interpretation 
in terms of the dispersion relation for a mode propagating through a medium is possible 
only if the imaginary part of uj k is small compared to its real part. In this case the mode 
can be visualized as a particle (or antiparticle) with an energy u r and a damping 7. Under 
such assumptions, each one of Eqs. (|2.8| ) and ( p.9|) yields two distinct solutions, one with 
positive energy and the other with negative energy, whose physical interpretation has been 
discussed in detail in Ref ||16|| . Here we will concentrate on the solution of Eq. (|2.8|) having 
a positive real part, which corresponds to the neutrino mode with energy u r , but similar 
considerations and results apply to the other solutions. 

It is convenient to make the decompositions a = a r + Wj and b = b r + ibi. Then, using 



Eqs. ( |2.10| ) and (|2.11|) , expanding in powers of 7 and retaining only terms that are at most 
linear in 7, aj and bi, from Eq. (|2.8|) we obtain 



f r (u r ,K) = 0. 



(2.12) 



and 



7 
2 



(2.13) 



with 



f r = (1 - a r ){uj - k) - b r 
fi = -ai{u — k) — h. 



(2.14) 
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Only approximate analytical solutions of Eq. ([2.12 ) are known |TJ|. At the one- loop level 



both b r (ui r , k) and a r (u r , k) are of order g 2 /M^. Therefore, to this order the energy of a 
massless neutrino is 

u r = K + 6 r (^ r ,K)| Wr=K . (2.15) 

On the other hand, as we will show, the imaginary part of £ vanish to 0(g 2 ) and to this 
order there are not contributions to lmu K . For the perturbative solution of Eq. ([2.1 3|) 
around w r = k, we have 

l^- bi (co r ,K)l r=K , (2-16) 

with bi{u r ,K) of 0(g 4 /M^). 

The matter effects on the neutrino oscillations are conveniently incorporated by means 
of the effective potential V. This can be introduced by subtracting the (vacuum) kinetic 
energy from the real part of the neutrino dispersion relation ||: 

V = u r -K = b r (uj r ,K)l r=K . (2.17) 

In the literature it is also customary to use a refraction index, which is defined by 

n(«) = — . (2.18) 

In the approximation we are working on, and utilizing Eqs. (|2.11| - |2T7|) , it follows that its 
real and imaginary parts are related to the effective potential and the damping rate by 

V 

Ren(K) = 1 , 

K 

lmn(n) = —, (2.19) 
2/x 



with V and 7 given by Equations ( |2.17| ) and ( p,16| ), respectively. 



III. EFFECTIVE POTENTIAL 



According to Eq. ( |2.17| ), for a perturbative solutions of the dispersion relation around 



the vacuum value uj r = k, the effective potential is given by the real part of the coefficient 
of i/l in the neutrino self-energy. This is only true in the lowest order, in general, V will also 
receive contributions from a r . From Eq. (|2.5| ), it is easy to see that 



b r (u,K) = ^Tr{/ReS} . (3.1) 



with 

1 



X = —(u U - kM) . (3.2) 

K 



f5 



In Ref. H the real part of the neutrino self-energy was calculated in a general gauge up 
to terms of order g 2 /M^. It was shown, that although the self-energy depends on the gauge 
parameter, the dispersion relation is independent of it. Taking this result into account, for 
simplicity we will work in the unitary gauge. Furthermore, for physical situations where the 
temperature is much lower than the masses of the gauge bosons, the propagators for the W 
and Z can be taken the same as in the vacuum, namely, 

with B = W, Z. 

We shall assume that neutrinos are massless, so at the one-loop level the only contribu- 
tions to £ arise from the diagrams depicted in Fig. 1. 

For a neutrino v\ (I = e, fi, r) that propagates through a medium with a momentum k, 
we split the different contributions to b l according to the processes in Fig. 1 as follows: 

K = b\ ad + b l z + b l w , (3.4) 

corresponding to the tadpole, Z-exchange, and W-exchange contributions. In this case, the 
background dependent parts of each term in the right hand side of the previous equation 
can be worked out as 



jTr{^ ia }Af (0)X; / ^r 3 ^{^ + mf)r%}5( P 2 -m}) Vf (p.u), 



.J4 

°tad 



(27T) 

b l z = ~\j ^Tr {Jl^r*,} Af(k - p)5(p 2 ) rj^p ■ u) , 

b l w = ~\j ^Tr + mi )Tj} A#(* - p)5(p 2 - m]) m (p ■ u) . (3.5) 

In by/ the charged lepton in the internal line has the same flavor that v\ , while in the 
tadpole contribution btad the summation is over all the fermion species present in the thermal 
background. In the previous expressions the vertices are given by 

With the vector X\ and axial couplings Y\ given for the charged leptons by 

X t = ~\ + 2 sin 2 9 W , 

Yi = \\ (3.7) 



for the neutrinos by 

1 
2 



and for the nucleons by 

1 

2 



X p = ^ — 2 sin 2 



y p = X n = -Y n = -I . (3.9) 



According to Eq. ( p. 171 ) b l r has to be evaluated at u = k, in this case the quantity k^. in 



( ft.2|) reduces to the energy projector for a massless particle in the vacuum and may be 
replaced by its usual expression in terms of the free spinors. In general for the external lines 



the spinors to be used are the solutions of the Dirac equation in the medium [13]. However, 
within the approximation we are using, they can be approximated by the vacuum solutions. 
For the internal fermionic lines we substitute j> + m/ by its corresponding energy projector. 
It is now straightforward to show that any of the b l in Eq. Q3.5|) can be written in the the 
form 

b l r (u,K)l= K =±;(M l f ), (3.10) 

where M. 1 ^ is the tree-level invariant amplitude for the forward scattering v\j — > uif . The 
corresponding Feynman diagrams are shown in Fig. 2. Notice that these diagrams are 
obtained from the self-energy diagrams in Fig. 1, if we cut along the internal fermion line. 
Brackets in the previous expression represent the thermal average given by 

<X/> = E/ ^Stf-™ a f)M l f T, f (p). (3.11) 



The previous result shows, that utilizing the methods of finite -temperature field theory (at 
the g 2 order), the neutrino effective potential reduces to the thermal average of neutrino 
forward scattering amplitude. In fact, it is interesting to combine Eqs. ( |2.17 ), (|2.19|) , and 
( p.lO| ) to write the real part of the index of refraction as 



2tt 



n = l + — (M l f ) . (3.12) 



That generalize the zero temperature result |17| n = 1 + 2tt M/k 2 , by simply adding the 



thermal average of forward scattering amplitude, a result that proves that background does 
not spoil the coherent condition of the forward scattering processes. 

In the rest of this section we derive the effective potential for a neutrino propagating 
in a thermal background composed of electrons, protons, neutrons, neutrinos, and their 
respective anti-particles. As previously mentioned, the Feynman diagrams in Fig. 2 are 
obtained cutting the self-energy diagrams of Fig. 1 along the internal fermion line. Diagram 
(a) is obtained from the tadpole self-energy, the result is the same for any neutrino flavor, 
and one has to sum over all the fermions / present in the background. Diagrams (b) and 
(c) correspond to cutting the Z- and W-exchange self energy diagrams, consequently the 
background fermion necessarily has the same flavor as the test neutrino. In this way, AiK 
can be write as 

M l f = M a + M b 5 fl/l +M c 5 fl , (3.13) 
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with 8u = 8 vm = 1, 8ft = for / ^ Z, and 8 fui =0 for / ^ v x . Here 



M a = -ju^k, s')Tf a u vl (k, s')Af(0)u f (p, s)TfpU f (p, s) , 
M b = ^u Vl (k,s)T z fa u f (p, s)Af(k - p)u f (p, s)T z fl3 u{k, s') , 
M c = ~u n (k, s')T^u f (p, s)A$(k - p)u f (p, s)Tju Ul (k, s') 



(3.14) 



Since we are interested in contributions to V of order g 2 /Myj, we expand the gauge 
propagator in power of M^ 2 up to the second order 



9af3 



9a/3 



M 2 B 



q 2 -M 2 B M 2K 
Using this expansion and neglecting quantities of order m 2 /Mg, we find 

M a ^2V2G F X f k-p, 



(3.15) 



M b ^V2G F [k-p 
M c « y/2G F [k ■ p 



2(p ■ kf 

2{p ■ kf 

M 2 W 



(3.16) 



where q = p — k and Gpj\f2 = g 2 /8Myy 2 is the Fermi coupling constant, and the factors 
Xj are given in Eqs. ( |3.7| ), ( |3.8| ), and ( |3.9| ) for leptons, neutrinos and nucleons respectively. 

Once the previous results are substituted in Eq. ( p.ll|) , the different contributions to V 
can be expressed in terms of the following integrals 



d A p 



p fl 8(p — m )rjf(p ■ u) = A 



d A p 



p»p"8(p 2 - m 2 )7] f (p -u) = B uV + Cg** v . 



(3.17) 



The scalar quantities A, B, and C are easily evaluated in the rest frame of the medium; the 
results are 



.4 
B 

C 



m}((±-)N f + (i-)iV>) - ((E f )N f + (E f )N^ 
m){^)N f + (i-)JV>) - 4((E f )N f + (E f )N^ 



(3.18) 



In these equations Nf (Nf) represents the density of fermions (anti-fermions) in the back- 
ground 



N 



fj ~ 9f 



d 3 p 



(2tt) 3 e^fT^f) + l 



(3.19) 
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and (Elf) ( j — 1> — 1) denotes the statistical averaged of £?/ and 

Here, -E/ = y'p 2 + m'j and (7/ is the number of spin degrees of freedom (g u = 1 for chiral neu- 
trinos and gf = 2 for the electron and nucleons). Collecting these results it straightforward 
to write down the effective potential for a neutrino of a given flavor. 

The contributions to the effective potential for the various neutrino flavors and back- 
ground particles are listed in Table I. They agree with the results given in the literature 



Neutrino 


Background Particle 


Veff 




P 


±^(1- Asm 2 9 w )(N p -N p ) 




n 




V e 


e 


±^§(1+4 sin 2 9 w )(N e -Ne) 
+ 2 ^u£ ml m±)+N e (±)-^(N e (E e )+N e (E e ))] 






±^(^ e - ~ + N, e {E 9e )] 




e 


±^§(-1+4 sin 2 9 w )(N e -Ne) 









Table I. Effective potential for a neutrino propagating through a medium. The + (— ) sign refers 
to neutrinos (anti- neutrinos). 



IV. IMAGINARY PART 



The discontinuity, of the self energy is related to the damping rate 7 that determines 
the imaginary part of the dispersion relation (see. Eq. 2 . 1 If ) , additionally the damping 



rate can be interpreted as the rate at which the the single-particle distribution function 
approaches the equilibrium form [|ll]]. The former interpretation follows if one consider a 
particle distribution that is slightly out of equilibrium, hence one has 

^ = -/r a + (i + a/)r c , (4.1) 

where T a and T c are the absorption and creation rates of the given particle respectively. The 
parameter a distinguish bosons (a = 1) and fermions (er = —1). The previous equation has 
for general solution 

/ K, t) = Tc + C K) e^-^ , (4.2) 

where C (u r ) is an arbitrary function that dose not depend on time. Creation and absorption 
rates are related by the KMS relation 
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T a (uj) = e^ T T c (uj) 



(4.3) 



Consequently Eq. ( |4.2|) can be written as 



/ (w r , t) = — ^ + C ( Wr ) e" 2 ^ , (4.4) 

where 7 = (r a — oT c )/2 is defined as the damping rate. Therefore, 7 can be interpreted as 
the inverse time scale it takes for a thermal distribution to reach equilibrium. The sign of 
the damping rate must necessarily be positive for an stable systems. Additionally, the form 
of the dispersion relation implies for a normal mode to propagate that 7 is small compared 
to u r . For neutrinos this conditions is satisfied in normal matter, such as the core of the 
Sun. However, in a CP-symmetric medium the leading contributions to the real part vanish. 
Thus, the first non- vanishing contributions are of order Gp/M^; under this circumstances 
7 can become of the same order. 

To obtain the neutrino damping rate we have to evaluate S 12 and then use Eqs. ( |2.2j ),(f2~5| ) 
and (|2.16|) . As explained further ahead the one loop contributions to £ 12 cancel. The 
diagrams that contribute to S 12 at the two loops level, are depicted in Fig. 3; there are also 
diagrams similar to those on (c) with the W and Z lines interchanged. According to the 
Feynman rules on the real time formulation of the FTF, the contributions of diagrams (a) 
and (b) can be written as 

d A q d A p 



i^(k) = -J (2 ^ 4f2 ^ 4 r^(^ 12 (g))r-Tr{r?(^ 12 ( P )) 



(2tt) 4 (2tt 

" MU " \„,l" - '/IHI'M,,; 

Similarly, for diagrams (c) and (d) we have 



xrf (iS 21 (p + q - k)MA*(k - g))n*(A* (* - q)) 22 . (4.5) 



-ffii a (*) = / (2^)4 ^y T i^M)^s 21 { q + P - fc))r?(»g 12 (p))rs 

xi(A* a (k - q)) u i(A* p (k - p)) 22 . (4.6) 

In the above expressions A, B = W, Z and = — T", with T° representing any of the vertex 
or rfj given in the previous section. The Su and 5* 2 i components of the propagator 
matrix of the fermion are given by 

Si 2 (p) = 2ni5(p 2 - m 2 )[r]f(p) - 9(-p ■ u)}(-/> + m) , 

S 2 i(p) = 2ni5(p 2 -m 2 )[r] f (p) - 6(p ■ u)](j) + m) . (4.7) 



In principle, the internal vertices should be added over the thermal indices (a = 1,2). 
However, since temperature is small as compared to the gauge boson masses, the matrices 
of the W and Z bosons propagators are diagonal with (A^g(/c)) n = — (A^/c))^ = Avg, 
where Afg is the vacuum propagator given in (|3.3| ). This explains the cancellation of the 
g 2 contribution to the neutrino damping rate. The one loop contribution to E12 is given 
by diagrams similar to those in Fig. 1 with iS\ 2 and i(A^ /3 (A;))i 2 replacing the internal 
fermionic and bosonic lines; however the bosonic propagator is diagonal, hence S 12 cancel 
at this order. 

in 



According to Eq. ( |2.16| ) 7 is directly proportional to bi(u, k), that is given by 



bi(u,K) = ^Tr{^ImS} , (4.8) 

with ft defined in Eq. (|3~2|) . 

Taking into account the previous results it is demonstrated after a lengthy calculation 
that the neutrino damping rate can be expressed in the form 

^JJ^Lww + c zz _ y d ab\ (49) 
where C AA is obtained from Eq. ( |4.5| ): 

cAA 4 / (§ ^§)^ m{i + + "*> 

x (2tt) 4 5 (4) (Q -(q + p- k))5(q 2 - m\)5{p 2 - m 2 2 )5(Q 2 - mj) 
x [VF(q x ) ~ 0(-q ■ u)] [v F (p y ) ~ 0(-p ■ u)} [r] F (Q z ) - 6(Q ■ u)] , 

(4.10) 

and D AB is obtained from Eq.( [4.6D : 

^ = \ I ^^^{^™ + + 

xTO + m 3 )Tl }i(A A a (k - g)) u i(A5,(fc - p)) 22 ] 

x (2tt) V 4) (Q -(q + p- k))5(q 2 - m\)5(p 2 - mj)5(Q 2 - m 2 3 ) 
xfoHfe) - 0(-q ■ u))[rj F (p y ) - 6(-p ■ u)} [r) F (Q z ) - 9(Q ■ u)) . 

(4.11) 

For later convenience an extra integration over the momentum Q has been introduced. In 
what follows we shall see that the neutrino damping rate can be expressed in term of ampli- 
tudes squared and weighted with the statistical factors that account for the various physical 
processes. To derive these results we notice first that fermion propagators in equations 



and (|4.6| ) are either type 12 or 21. According to Eq. ( |4.7|) the propagator Su(p) contain a 
delta function S(p 2 — m 2 ) and a factor (p 1 + m). The delta function put the fermion on their 
mass shell mass, in other words self-energy diagrams in Fig. 3 are cut along all the internal 
fermion lines. Whereas the second factor is concerned, insertion of the fermion projectors 

p 1 + m = U (.P, s )u(p, s) 

s 

p 1 — m = v(p, s)v(p, s) , (4-12) 
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allow us to rewrite the resulting expressions in terms of amplitudes for the physical processes 
arising from the cuts. The bosonic lines are not cut because they do not include thermal 
distributions for T <C M w . 

For definitiveness let us consider diagram (b) in Fig. 3, and also that the fermion in the 
internal loop is a proton (/ = P). When the diagram is cutted as shown in the figure, we 
obtain a series of physical processes for the neutrinos vi, v-i and protons P\ and Pi. Of these 
particles one of the neutrinos (y*) is considered a test particle, all the other are thermalized. 
According to the notation in Eq. [4.6| , the momentum and chemical potentials are assigned 
as: (ui : k,fi vi ), [y 2 : q,fi U2 ), (Pi : Q, fM Pl ), and(P 2 : P,Hp 2 )- With momentum and charge 
conservation conserved depending of the process, e.g for V\P\ — > v^Pi. 

k + Q = q + p, 

Pvt. + At = Hvx + lip x ■ (4-13) 

The processes obtained from the mentioned cut rules include the two neutrinos and the two 
protons distributed into the initial and final states in all possible ways. Hence, in general we 
expect to obtain 16 different processes; this is explicitly display in Eq. ( [4.16 ) The resulting 



expression come out with the appropriated thermal distribution, for this we have to rewrite 
the thermal contributions that appear in Eq. ( |4.10 ) utilizing the following identities: 



T]f(xk) - 0(±k ■ u) = =Fe(A; • u)n f (^x k ) , 
n f {x k ) = e~ Xk n f (-x k ) = e~ Xk [l - n f (x h )) , 

n f (x k )n f (x q )n f (x p ) = n f (x k )n f (x q )n f (x p ) 



n f (x Q ) 



+ e-*>—+*> (1 - n f (x k )) (1 - n f (x q )) n f (x p ) , (4.14) 



where: 



x k = (3{k ■ u - n Vl ) , x q = (3{q ■ u - fi V2 ) , 

Xp = @(jp • U - fJL Pl ) , Xq = P(Q ■ U - (JLpz) . (4.15) 

Taking into account these consideration and performing an integration over the time-like 
components of the momentum integration, it is possible to cast the contribution to 7 arising 
from the proton loop in diagram (b) of Fig 3 into the following form 



k I ^hyiSky w^t' 4 ^' +Q - 9 ~ p) 1 M{Pv " Pv) ? 

x [(1 - n v {E q )){l - n P {E p ))n P {E Q ) + n u {E q )n P {E p ){l - n P {E Q ))} 

+5 {A) (k + Q + q-p) \M{PW^P) | 2 [n v {E q )n P {E Q ){\ - n P {E p )) 

+ (1 - n p {E q )){l - n P {E Q ))n P {E p )} + 5 {4 \k + Q + p - q) \ M{PPu <-> v) | 2 

x [(1 - n v {E q ))n P {E p )n P {E Q ) + n v {E q ){\ - n P (E p ))(l - n P {E Q ))} 

+5^(k + Q+p + q) I M(PPvv^O) | 2 [n p (E q )n P (E p )n P (E Q ) 

+ (1 - n s (E q ))(l - n P {E p )(l - n P {E Q ))] + 6&{k -Q-q-p)\ M(u <-> PPv) 

x [(1 - n v {E q )){\ - n P {E p )){\ - n P {E Q )) + n v {E q )n P {E p )n P {E Q )} 
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+6^(k + q -Q- p )±\M(W^PP)\ 2 [n P (E q )(l - n P (E p ))(l - n P {E Q )) 

+ (1 - n D {E q ))n P {E p )n P {E Q )} + 5 (4) (k+p — Q - q) \ M[Pv <-> Pv) | 2 
x [(1 - n v {E q ))n P {E p ){\ - n P (E Q )) + n u (E q )(l - n P (E p ))n P (E Q )] 
+8^(k+p + q-Q) | M(PuV^P) | 2 [n v {E q )n P (E p ){l - n P {E Q )) 

+ (l-n D (E q ))(l-n P (E p ))n P (E Q )} 



(4.16) 



Regardless of its length the interpretation of this equation is quite simple. The first two terms 
are interpreted as the absorption and emission of a neutrino via the dispersion vP — > vP 
with statistical weight (1 — n u ){l — n P )n P and n u n P {l — n P ), respectively. As expected a 
rif factor appears for each background fermion in the initial state, whereas fermion in the 
final state contribute with a Pauli blocking term 1 — rif. As already discussed for femions 
the absorption and emission decay rates must be added |]TT|. The scattering amplitude for 
both processes are the same, because they are related by CPT inversion. Similarly the third 
and fourth contributions represents neutrino annihilation and creation via the PvV — > P 
and P — > PvV respectively; they include, as expected, the statistical factors n v n P {\ — n P ) 
and n P (l — n^l — n P ). The same reasoning applies to the remaining terms. 

Taking into account the 5— function constrains some of the quoted processes are not 
allowed. In what follows we focus on conditions with temperatures rrif « T <C My/ where 
for example the composition of the primeval plasma is dominated by (anti-) neutrinos (anti- 
electrons, nucleons and photons. We recall, see discussion below Eq. [4.4] , that for fermions 
the contributions to 7 of decay and absorption add together. Hence the statistical factors 
appearing in the previous equation can be simplified. For example the absorption v*P — > vP 
and decay vP — > u*P, where {y* is the test particle) add according to 

(1 - n u {E q )){\ - n P (E p ))n P {E Q ) + n v (E q )n P (E p )(l - n P {E Q ))} = (4.17) 

n P (E p )(l - n P (E Q ))} + n u (E q ) (n P (E p ) - n P (E Q )) . 

However we can drop out the last term in the right hand side of the equation, because its 
contribution vanish when substituted into Eq. ( |4.16| ). 

With all these results we finally find that the neutrino damping rate can be expressed as 

7 = ± M!(I/^ tl ^|A^e~,e)| 2 

1 2u {2 J (2^)32^(2^)32^ (2vr) 3 2£ Q 1 1 )l 



xd 2 (k + p-q-Q)(l- n e (E Q ))n e (E. 



vi 



+ all the possible processes j , (4-18) 

where the + (— ) sign stands for test neutrinos (anti- neutrinos). For all possible processes we 
mean all the kinematically allowed processes obtained by cutting all the fermionic internal 
lines of Feynman diagrams shown in Fig. 3. The corresponding processes and their cross 
sections are listed below in (|4.22|) . 



The damping rate can be written in terms of the thermal average of the cross section 
For the dispersion vf «-> uf, the differential cross section is given by 
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where VJ. e i is the relative velocity between the neutrino and the background fermion; as we 
are considering massless neutrinos we simply have V re \ = 1. 
Thus, 7 reads 



7 



±£«//^<^/> ( 4 - 2 °) 



where a/ = 1/2 for f = e,n,p and a/ = 1/4 for / = v and 

<**> = M^J mA, ^ .-«(!- »/(«,))■ (4.2D 

is the cross section thermally averaged by the Pauli blocking term (1— rif(EQ)). In Eq. (|4.18 ) 
the ellipsis represent the other possible process, in each case the corresponding statistical 
factors in Eqs. ( 4.18 , |4.19| ) and the dispersion amplitude are replaced by the pertinent 
factors. We recall that according to Eq. ( p.!9| ) 7 is directly related to the imaginary part 
of the index of refraction. Hence, Eq. (|4.20|) can be identified with the optical theorem. 

In what follows we shall apply the previously obtained results to explicitly evaluate 
the neutrino damping rate in a background composed of neutrinos, electrons, protons and 
neutrons. We suppose that the Pauli blocking term can be neglected, in addition we consider 
temperatures nif ~ T <C Mw, hence in the thermal averages we can assume that q 2 -C M^r. 
First let us consider the cross section for the relevant processes. It is common to quote 
the cross sections, assuming ultrarelativistics neutrinos and neglecting the fermion masses. 
However for conditions as those of the early universe, temperature and consequently the 
average neutrino energy can be comparable to the nucleon masses, and sometimes to the 
lepton masses. Hence, keeping fermion masses, the various neutrino cross sections can be 
calculated as 

771 ^ S 777 ^ 

a(u e e <-> u e e) = 165 2 + 125 + 3 - (405 2 + 305 + 6) + (125 2 + 85 + 1) e 



s 2 



(165 2 + 26)^$ + 45 2 ^| 
s 6 s 4 



2 771 ^ ^ 777 ^ 

a{v^ T e <-> i/ MT e) = 165 2 - 125 + 3 - (265 2 -155 + 3) — e - + (125 2 - 85 + 1)- e 



s 2 



( 85 2_ 35) ^| + 45 2^e 



S 3 S 4 



777 ^ *^777 

a{y e e <-> u e e) = 165 2 + 45 + 1 - (405 2 + 105 + 1)-^ + (125 2 - 45) 



s s 2 

- (165 2 + 105 + 1)^ + (25 + l) 2 ^ , 
s 6 s 4 

777^ StT?^ 1 

a(u„ T e <-> u u T e) = 165 2 - 45 + 1 - (405 2 - 165 + 1)— ^ + (125 2 - 45)^ 

s s 1 



- (165 2 - 125 + 1)^ + (-25 + 1 



s 3 v y s 4 
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diviVi <-> ViVi) = 8, 

a(uiUj <-> UiUj) = 6, 

biviVj <-> z/jZ/j) = 2, 

-/ \ 6m^ 3m 
aiVin <-> vm) — o h 



s s 2 
s s 3 




2 4 
777 771 

<-> M = (165 2 - 125 + 1) - (405 2 - 305 + 6)— ^ + (365 2 - 245 + l)-£ 

s s 2 

_(165 2 -65)-f + (45 2 -25)^, 

s 6 s 4 

2 4 
777 777 

a(uip <-> i/ip) = (165 2 - 45 + 1) - (405 2 - 105 + 1)— ^ + (365 2 - 12<f)-£ 

s s 2 

777 u T77 

_(l 65 2_ 105 + 1) !^ +(45 2_ 25) ^ ) 
S d S 4 

/ 4777 2 \ 1/2 / 1 Am 2 " 

~o{v e v e ^ ee) = ( 1 - -p J [ 165 2 + 85 + 2 + (85 2 + 45 - -) ' 

<- ee) = fl - ^ 7 ( 165 4 - 85 + 2 + (165 2 - 85 + ^)=^ ) , (4.22) 

where a = a/ao with o"o = j^s, i = e, fi,r, s = (k + p) is the Mandelstam variable, and 
5 = sin 2 9w ~ 0.229. These results reduce, in the zero fermion limit, to those found in 
Enqvist, Kainulainen and Thomson || and Langacker and Liu [p~C}1 

Once the cross sections are inserted in Eq. ( f4.20|) the thermal averages should be evalu- 
ated according to the constrains of the problem. If we consider temperatures well bellow the 
nucleon mass, then the proton and neutron contributions will be suppressed by the Boltz- 
mann factor, and their contribution neglected. On the other hand if T pa rrif the complete 
average of the cross sections in Eq. (|4.22|) should be considered. In what follows we consider 
situation in which rrif < T and we retain terms of order m'j/s in ([4.22 ). This lead us to 
consider thermal averages that contain integrals of the following type 

^3 (* • p) (k ■ u) 9{k ■ u)5(p 2 - m))n f {p) . (4.23) 

These integrals are easily evaluated in terms of the thermal average of the fermion density 
and energy, utilizing the Eqs. ( |3.17| - |3T20| ). The results can be collected in a general formula 
that gives the contributions to the neutrino damping rate arising form various background 
particles (to first order in rrif): 



£ {2uA f (E f } + m) (A f + Bj)} , (4.24) 



12tt Y g f 



where u is the neutrino energy, summation in / is taken over fermions in the background, 
and the corresponding factors Af and Bf are summarize in Table 2. In this equation Nf 



15 



represents the density of fermion of antifermions in the medium (Eq. [3.19| ) and the statistical 
averages of Ef is defined in Eq. ( p.20[ ). 



iieuu mo 


rjacKgiomiQ ienmon j 
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10 -f- LZO "T O 


('aha 2 _j_ ^nA _l f\\ 
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1 fi A 2 1 9 A 4_ 3 


^9A 2 ^flA 4- 
— ^ozo — OUO T Ul 


Ve 


e~ 


165 2 + 45 + 1 


-(405 2 + 105 + 1) 




e~ 


165 2 -45 + 1 


-(405 2 - 165 + 1) 


Vi 


Vi 


12 





Vi 


vj, i ^ j 


6 





Vi 


Vi 


8 





Vi 


Vj, i ^ j 


2 





Vi 


n 


3 


-6 


Vi 


n 


2 


-2 


Vi 


P 


165 2 - 125 + 1 


-(405 2 - 305 + 6 


Vi 


P 


165 2 - 45 + 1 


-(405 2 - 105 + 1) 


z/ e 77 e <-> ee 


annihilation process 


165 2 + 85 + 2 


-6 


v^rV^r <-> ee 


annihilation process 


165 2 - 85 + 2 


-1 



Table II. Coefficient Af y Bf appearing in Eq. ( gjg) for various processes between neutrinos and 
the quoted background particles. Here 5 = sen 2 9\y. 



The next step is to quote some explicit results. Consider a CP-symmetric plasma com- 
posed the three type of neutrinos v e , and u T , electrons and their corresponding antipar- 
ticles. Considering test particle, the contribution to the neutrino damping rate 

arising from the background neutrinos is given by 

G 2 , 



7 



7T" 



(4.25) 



where u is the neutrino energy. Whereas the contribution of the electron and positron 
background to the v e damping rate yields 



G 2 

7 = 0.39^ wT 4 /) 

7T 



T 



T 



-0.6- 



(4.26) 



where the functions p and e are defined by 



p(0 

c(0 



cfx — . 

ev/^+i 

oo 



dx- 



X 



(4.27) 



p(0 e ^x 2 H 2 + i 

Similar expressions can be obtained for the other processes. These results reduce to those 
in references [||,||[10| in the limit of zero electron mass. 
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V. CONCLUSIONS 



In this paper we have systematically derived the neutrino damping rate in real-time 
thermal field theory. Starting from the discontinuity of the neutrino self energy at the two 
loop level, we prove that the damping rate can be expressed as integrals over space phase 
of total cross sections, weighted with statistical factors that account for the possibility of 
particle absorption or emission from the medium. 

Cutkosky rules are used to obtain the self-energy imaginary part at zero temperature. 
Weldon 0] and Kobes and Semenoff [|l^] (see also have studied the corresponding 
Cutkosky rules at finite temperature. In these references it is shown that for certain specific 
examples at one loop order the discontinuity of the self-energy can be expressed in term of 
amplitudes squared and weighted with the statistical factors that account for the various 
physical processes. Here we prove that these result stand valid for the neutrino damping 
rate at the two loop order. 

The complete results that account for all possible processes that contribute to 7 appear 
in Eq. ( 4.16| ). Depending of the physical conditions some of these processes are forbidden. 



Specific results for conditions such as those of the early universe where the primeval plasma 
is composed of (anti-)neutrinos, (anti-)electrons, and nucleons were obtained. For those 
conditions the fermion masses are not always negligible, consequently we report a general 
formula (Eq. 4.24]) that includes mass correction to first order in mj/(s); however, further 



improvements are easily obtained utilizing the values for the cross sections in Eq. ( 4.22 ). Our 



results, summarized in Eq. ([4.24 ) and Table II, should be useful for the study of neutrino 



processes in the early universe, as well in some astrophysical scenarios. 
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FIGURES 



C)f z w 




FIG. 1. One- loop diagrams for the self-energy of a neutrino in a thermal background of charged 
leptons nucleons and neutrinos. In (a), / stands for any fermion in the background. In (c), the 
charged lepton I is of the same flavor as the neutrino. 
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u t v l v l v l v l I 

(a) (b) ( C ) 

FIG. 2. Tree-level diagrams for v\j — > v\j forward scattering; / represents the background 
fermions. These diagrams are obtained when the corresponding self-energy diagrams [Fig. 1] are 
cutted along the internal fermion line. 
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(c) z (d) z 

FIG. 3. Two-loops contributions to the self-energy of a neutrino in a thermal background of 
charged leptons, nucleons, and neutrinos. 
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